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Cross-Reference 

This patent application claims the benefit of United States Provisional Patent Application 
Number 60/050389 filed on June 20, 1997. 



This invention relates to the general subject of seismic exploration and, in particular, to 
methods and devices to measure seismic anisotropy in-situ, with improved accuracy and resolution. 



It is well-known to those familiar with the art of exploration seismology that the rocks of the earth's 
subsurface exhibit seismic anisotropy, and that failure to account for this can confuse the interpretation of 
seismic data. For example, surface-reflection seismic data acquired by experiments in anisotropic rock, 
but interpreted via isotropic theory, can yield incorrect estimates of reflection depths, incorrect analysis of 
Amplitude-vs.-Offset (AVO) measurements, etc. 

Seismic anisotropy is that property of materials by which the velocity of seismic waves (taken 
here as P-waves) depends upon the direction of travel of such waves. Polar anisotropy is that special 
case where the variation is primarily or exclusively with the poiar angle (i.e., the angle from the vertical, 
rather than the azimuthal angle). It can be caused, for example by preferred alignment of rock grains {as 
in shales), or by sedimentary layering on a scale much smaller than the seismic wavelength. Because 
such layering is common in the earth's crust, and because seismic waves of differing wavelength are 
affected differently by such layering, the anisotropy measured by any seismic experiment depends upon 
the wavelengths used in the experiment. Hence, experimental results from widely differing wavelength 
bands may not be combined to estimate the anisotropy (more on this later). Since subsurface rocks in 
different places possess seismic anisotropy to varying degrees, accurate local measurements of the 
anisotropy are routinely required in order to make corrections, thereby avoiding errors such as those 
mentioned above. 

FIG. 1 is a good starting place to illustrate the problem. In a typical seismic experiment, seismic 
waves are emitted by a source at S, reflected by a rock-layer at depth, z ,and received by a receiver at R, 
separated from S by a distance x. With the geometry shown, and if the rock is uniform and isotropic, it is 
well-known that the arrival time t(x) of the wave arriving at offset x is given by: 
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t 2 (x)=t 2 0 + 
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where to is the vertical travel time or the arrival time (of the wave) at the special place where x = 0 (i.e., the 
place of normal incidence), and where V is the velocity of seismic waves. The quantities V and to are 
conventionally determined by measurements of t(x), with V determined as the parameter which best fits 
the hyperbolic variation of arrival time with offset. The (initially unknown) depth z to the reflector is given 
by: 




(2) 



With V and to measured, the depth z is thereby determined. 

Of course, in any practical situation, the geometry will be more complicated than that shown, and 
the material properties will not be uniform. For example, FIG. 2 illustrates the situation where the reflector 
is dipping at angle (j>. It is well known (See Slotkin, M. M., "Lessons in Seismic Computing," Soc. Expl. 
Geoph., 1959, Tulsa, Oklahoma, the teachings of which are incorporated herein by reference) that: 



where to is given by Equation (2). However, the simple situation shown in FIG. 1 suffices to illustrate the 
difficulties which arise if the rock is considered to be anisotropic. 

If the rock is anisotropic, then Equation (1) no longer applies as written. To generalize it for 
anisotropy, one has to specify what sort of anisotropy the rocks possess. In the simplest anisotropic case, 
the velocity varies only with the incidence angle 0 shown in FIG. 1, and not with azimuthal angle. In this 
(i.e., polar anisotropic) case, it is well-known that the angular variation of velocity depends upon four 
material parameters, following a complicated formula known to those skilled in the art. 

However, it is also well-known that, if the anisotropy is "weak," then the P-wave velocity V p (6) 
depends upon only three parameters which may be chosen (modifying Thomsen, "Weak elastic 
anisotropy," Geophysics . 1986, 51(10), pp. 1954 - 1966, the teachings of which are incorporated herein by 
reference) as: 



t 2 (x)=t^ + | r+ (^)sin«|, 



(3) 




(4) 
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10 



20 



where V 0 is the velocity in the special case that 0 = 0 (i.e., it is the vertical velocity), and 6 and r\ are two 
non-dimensional anisotropic parameters. Materials with "weak" polar anisotropy have both 5 and r\ much 

less than one. The horizontal velocity V p {9 = 90") is given by: 

V p (0 = 90°)= V 0 [1 + S+rj] (5) 

so that 

s = 5 + T| < 6 ) 
is the horizontal anisotropy parameter. The parameter 5 controls the near-vertical variation of velocity, 
and r\ is called the aneliiptic anisotropic parameter. 

With this notation, one can generalize Equation (1) to the Vertical Traverse Isotropic case: 



t 2 ( x ) =t o + vr- (7) 

mo 

for small offsets X < Z . The parameter V mo , called the "moveout velocity" or dx/dt, is given in 
15 terms of quantities previously defined by: 



25 



V mo = V 0 VT+25 = V 0 [l + 5] (8) 



Further, the unknown depth z is given, in terms of the vertical velocity V 0 or dz/dt, by: 



V„t 

z = 



o l o (9) 



(compare to Equation (2)). Those skilled in the art will recognize that the vertical velocity and the moveout 
velocity are features of the P-wave velocity function. 



However, the vertical velocity V 0 is nor given by the hyperbolic variation of t with x, which is V mo 
and which includes an anisotropic correction (see Equation (8)). Hence, although V mo is known, V 0 itself 
is not known, and so the measured vertical travel time to may not be interpreted to yield the depth. Thus, 
because of the presence of anisotropy, neither the anisotropy itself nor even the depth can be determined 
30 accurately. 
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If longer offsets, x > z , are acquired, then the arrival times are given by a formula more 
complicated than that given by Equation (7) {i.e., with an additional (non-hyperbolic) term involving the 
aneilyptic anisotropy parameter ri; See Tsvankin and Thomsen, "Nonhyperbolic Reflection Moveout in 
5 Anisotropic Media," Geophysics , 1994, 59(8), pp. 1290-1304, the teachings of which are incorporated 
herein by reference). This non-hyperbolic term has been specified (See Aikhalifah, "Anisotropy 
Processing in Vertically Inhomogeneous Media", Soc. Expl. Geoph. Expdd. Absts. , 1993, 65, 348 - 353, 
the teachings of which are incorporated herein by reference) in terms of T|, generalized as: 



10 However, this additional term helps determine the aneilyptic anisotropy parameter tj, rather than the 
vertical velocity V 0 and/or the near-surface anisotropy parameter 8, so the problem remains. 

Heretofore, the only practical seismic way to determine the vertical velocity V 0 has been to 
measure it (e.g., by drilling a hole, and by using Equation (8)). Thus, when both the vertical velocity V 0 

15 and the moveout velocity V mo are determined, the near-surface anisotropy parameter 8 may by found 
using Equation (8). This set of operations then determines all three parameters (V 0 8, r|) with seismic 
resolution (i.e., with the resolution of a seismic wavelength or more (e.g., several hundred meters)). For 
some purposes, this is sufficient, but for others it is not. For example, proper interpretation of the AVO 
effect requires an estimate of the near-surface anisotropy parameter 8 with much greater resolution. The 

20 AVO effect is often used to detect the presence of gas in the pore-space of the rock, but this detection can 
be confused by anisotropy. 

The argument is as follows: 

After certain "true relative amplitude" corrections have been performed, the amplitudes of the 
received signals in the reflection experiment of FIG. 1 may be interpreted in terms of the reflectivity. It is 
25 well-known (See Thomsen, "Weak Anisotropic Reflections, in Offset-Dependent Reflectivity - Theory and 
Practice of AVO Analysis," Castagna, ed., Soc. Expl. Geoph . . 1993, Tuisa, pp. 103 - 114, the teachings of 
which are incorporated herein by reference) that, if the contrast in elastic properties across the reflecting 
horizon is smalt, the reflectivity at small-to-moderate angles may be written in the isotropic case as: 





(10) 



30 



where 
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0 2Z 



(10a) 



and 



1 2 



( 



AV 
V 



f 2pO AG 



V 



G 



(10b) 



5 Here, Z is the average P-wave impedance across the horizon, V s is the average shear-wave velocity, 
and G is the average shear modulus. AZ, AG and AV are the corresponding jumps across the 

reflecting horizon. Anomalous values of R] result from the presence of gas on one side of the horizon, 
and this is a useful exploration signature. However, anomalous values may also result from anisotropy. 

The polar anisotropic version (See Thomsen, 1993, supra, and Rueger, "P-Wave Reflection 
10 Coefficients for Transversely Isotropic Media with Vertical and Horizontal Axis of Symmetry," Soc. Expl. 
Geoph. Expdd. Absts. 1995, 65, 278-281, the teachings of which are incorporated herein by reference) of 
Equation (10) is: 



15 



20 



R p (e') =1^ + 1?.; sin 2 e' 

where: 9' is the wave-front-normal angle; 



(11) 



R'=^<L ; and (Ha) 
2Z 0 



1 2 



AV, 
V 



2V S0 1 AG 



V v 0 ; G 



V v o v v o ' ^ 0 



£- + AS 



(11b) 



Here, Z 0 is the vertical P-tmpedance, V 0 is the vertical P-velocity, V s0 is the vertical shear velocity, G 0 
is the vertical shear modulus, and A signifies a jump at the reflecting horizon. R\ also contains the jump 
in anisotropy parameter 8 across the reflecting horizon. The ray angle 0' differs from the ray angle 9 by 
25 an anisotropic correction factor involving 5 (See Thomsen, 1986, supra ). 
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It is clear from Equation (11) that, if the anisotropic contribution is present but ignored, then an 
anomalous value of R[ , due to gas may be masked (i.e., a false negative) or mimicked (i.e., false 
positive) by the anisotropy. Therefore, in order to interpret J?( , in terms of gas, the anisotropic 
contribution must first be estimated and corrected for. However, this estimation must be locally 
5 determined, with a resolution of V* wavelength or less; hence, the process of Equations 8 and 9, offering 
resolution of more than one wavelength, is not suitable. 

Thus, a practical method and apparatus to do this with suitable resolution is needed. 

10 Summary Of The Invention 

Heretofore, no method was available to measure seismic anisotropy in-situ, with sufficient 
accuracy and resolution to solve many of the difficulties facing those skilled in the art. 

In accordance with the present invention an improverd method and apparatus is disclosed for 
measuring seismic polar anisotropy in-situ in a region of the earth that is characterized by geological 
15 formation elements and/or lithologic horizons. In one embodiment, the apparatus of the invention 
comprises: a recording system; a borehole having an axis that is deviated from the vertical by a known 
acute angle; a housing that is adapted to travel within the borehole and that is in electronic 
communication to the recording system, the housing carrying at least one source of acoustic energy and 
at least two receivers for receiving acoustic energy from geological formation elements and/or the 
20 lithologic horizons and from the source; and processing means for operating the source and the receivers, 
for recording the position of the housing in the borehole and for processing the recorded data in terms of 
both the direct raypaths from the source to the receivers and the indirect raypaths from the source through 
geological formation elements and/or lithologic horizons that are traversed by the source and the receivers 
so as to obtain measures of at least some of the seismic polar anisotropy parameters, the indirect ray- 
25 paths including reflections from geological formation elements and/or lithologic horizons that surround the 
borehole and that are at angles different from that of the axis of the borehole. 

In one embodiment, the processing means comprises means for measuring the anisotropic phase 
slowness over a two-dimensional suite of directions, and by constructing a slowness-surface plot, that is, 
a graph of vertical phase slowness versus horizontal phase slowness. This plot can also be interpreted as 
30 a plot of inverse phase velocity as a function of plane-wave propagation direction, by which one can 
measure phase slowness along the direction of the borehole. 

In another embodiment of the invention, a high resolution method is disclosed for measuring 
seismic polar anisotropy in-situ in a region of interest characterized by geological formation elements that 
reflect sonic energy. In that embodiment the process comprises the steps of: using a source to input 
35 sonic energy into the subsurface of the region of interest; recording sonic energy reflected from the 
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subsurface of the region of interest by using a plurality of sonic receivers that are located in a borehole 
that is located in the subsurface of the region of interest, wherein the borehole having an axis that is 
deviated from the vertical by a known angle; and interpreting the recorded sonic energy of the receivers in 
terms of both the direct raypaths from the source to the receivers and the indirect raypaths from the 
5 source to the receivers to obtain a measure of the seismic polar anisotropy parameters V 0 , T|, and 8. 

Numerous advantages and features of the present invention will become readily apparent from 
the following detailed description of the invention, the embodiments described therein, from the claims, 
and from the accompanying drawings. 

10 Brief Description of the Drawings 

FIG. 1 is an elementary diagram depicting a basic seismic reflection experiment; 
FIG. 2 is an elementary diagram depicting a seismic reflection experiment with a dipping horizon; 
FIG. 3 is a representation of conventional seismic reflection measurement tool in a vertical 
borehole; 

15 FIG. 4 is a representation of a seismic measurement tool located in a borehole which is not 

vertical and in which signals are processed in accordance with one aspect of the present invention; 

FIG. 5 is a representation of seismic measurement tool located in a borehole which is not vertical 
and in which signals are processed in accordance with another aspect of the invention; and 
FIG. 6 is a biock diagram of one embodiment of the invention. 

20 

Detailed Description 

While this invention is susceptible of embodiment in many different forms, there is shown in the 
drawings, and will herein be described in detail, two specific embodiments of the invention. It should be 
understood, however, that the present disclosure is to be considered an exemplification of the principles of 
25 the invention and is not intended to limit the invention to any specific embodiments or algorithms so 
described. 

The first step is to obtain seismic data in the form of seismic signal traces distributed over a three 
dimensional volume of the earth. Methods by which such data is obtained and reduced to digital form for 

30 processing as 3-D seismic data are well known to those skilled in the art (see FIG. 6). In one embodiment 
of the invention, a sonic logging instrument, sonde or tool 10 is used. The sonde 10 fits within a borehole 
9 and is in communication with a recording system 15 at the wellhead. The logging instrument or toot 10 
may be of conventional design, it may also be modified to suit the novel processing 14 employed by the 
present invention. Modifications may include a larger power supply (for greater range), a lower source 

35 frequency (for greater range), increased receiver spacing, additional receivers, tool centering devices, and 
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tool excentering devices. Referring to FIG. 3, a borehole source S of sonic P-waves transmits sonic 
waves which refract along the borehole wails 12 and are received at an array of receivers R-| through R n 
carried by the tool 1 0. 

The arrival times of the waves in a vertical borehole in a polar anisotropic formation are given by: 

t = t, + — (12) 

where t| is the arrival time at the first receiver R>| , and x is the offset measured from the first receiver R^ . 
The moveout is linear, rather than hyperbolic, as in Equation (7), since the wave travels in a straight line 
along the borehole walls 12 (i.e., no reflections). The vertical velocity V 0 is the average velocity over the 
array aperture (i.e., a few feet). In a conventional deployment (i.e., a vertical borehole), the waves travel 
in only one direction (i.e., along the axis of the borehole 9), so that no determination of anisotropy in-situ is 
possible. 

Referring to FIG. 4, in a polar anisotropic medium, if the borehole 9 is deviated from the vertical at 
angle 0, then the arrival times of the direct wave refracted along the borehole are: 



t=t.+ 7 — 7" (13) 

1 v p (e) 

where V p (9) is given by Equation (4). The present invention makes use of the other, indirect ray-paths 
shown in FIG. 4, using the same tool in the same context, to determine the anisotropy with higher 
resolution. 

In a deviated borehole, there will be reflected arrivals, as shown, coming from the ubiquitous 
horizontal layering. These arrivals have been utilized by others (Utard, M., C. Naville, and N. Rehbinder, 
"Study of Reflected Events Observed on Micro-seismograms Recorded With the EVA® Acoustic Logging 
System," 9th IntL Formation Evaluation Transactions . 1984, Soc. Prof. Well Log Analysts (Paris); and 
Fortin, J. P., N. Rehbinder, and P. Staron, "Reflection Imaging Around a Well With the EVA Full-waveform 
Tool", The Loo Analyst . 1991, May-June, pp. 271-278, the teachings of both are incorporated herein by 
reference) to construct "micro-seismograms" of the near-borehole environment. The "EVA" tool used by 
these workers was a logging instrument that was specially modified (e.g., in source-receiver offset, source 
characteristics, etc.) from more conventional designs, in ways that optimized the reception of these 
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reflections. However, velocity information, as in the present invention, was not used to estimate 
anisotropy. 

Those familiar with surface seismic reflections over dipping media will recognize that this 
geometry is closely related to that of the "dip-movement" problem (see FIG. 2, and Tsvankin, 1. "Normal 
Moveout From Dipping Reflectors in Anisotropic Media," Geophysics, 1995, 60(1), pp. 268-284, the 
teachings of which are incorporated herein by reference), with a similar result for the arrival times of the 
reflected waves: 



where the hyperbolic moveout velocity parameter V mo is given by Equation (8) (i.e., it contains the 



In accordance with one aspect of the present invention, data from the sonic tool 10 is treated 
exactly like surface seismic data to estimate this moveout parameter. For example, the data may be 
sorted into "Common Mid-Point (CMP) gathers," and maximum semblance methods may be used to 
determine the best value for V mo . 

As the sonde 10 moves through the borehole 9, it generates a set of source-receiver ray-paths 
that are similar in geometry to a moving single-ended shot spread in surface acquisition practice, and very 
much like a marine towed-streamer survey, if there are multiple reflectors nearby, it is even better than a 
towed-streamer survey. It is more like a group of towed-streamer surveys each done with different water 
depths, and even different water "tilts". 

The results may be processed directly using Equation (14), estimating the "Dip MoveOut" 
parameters V mo and t 0t via any of a number of conventional "DMO processes", with the angle 0 known, 
in this application, from the borehole trajectory. Also, the parameters V mo and t 0 may be determined by 
methods of pre-stack or post-stack migration, or reflection tomography. Either way, this yields a second 
independent equation for the three anisotropy parameters. 

A third equation is provided simply by Equation (9), with z known from the borehole trajectory. 




(14) 



anisotropic correction factor 
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These three Equations (9, 13 and 14) are usually sufficient to determine the three anisotropy 
parameters (V 0 , 5, and r\) with a resolution of a few meters, depending on the details of the geometry of 
the tool, the source characteristics, and the attenuation of the formation. 

Furthermore, if non-hyperbolic moveout from the reflections is detectable within the span of 
available source-receivers separations, then this non-hyperbolic moveout can furnish a fourth equation 
(See Tsvankin and Thomsen, 1994, supra ; and Alkhalifah, supra) that can be used to provide an over- 
determined set of four equations in three unknowns, which can be solved for a more robust determination 
of these unknowns, complete with error estimates. 

A More General Embodiment 

Thinking of a source and a receiver in homogeneous two dimensional space, one can imagine 
freely moving them about, and getting the direct arrival time for any combination of source and receiver 
positions. Since the medium is homogeneous, all that matters is the (vector) off-set from source to 
receiver; call that offset (Xq, Zq). 

For a given source-receiver offset (X 0 , Z 0 ), if one perturbs X 0 around that value, then one can 
directly measure the horizontal phase slowness for that (X 0 , Z 0 ). Similarly, by perturbing Z 0 , one can 
measure the vertical phase slowness for that (X 0 , Z 0 ). By doing that for several different values (X 0 , 
Z 0 ), one can directly construct a slowness-surface plot (i.e, a graph of vertical phase slowness versus 
horizontal phase slowness). This can aiso be viewed as a piot of inverse phase velocity as a function of 
plane-wave propagation direction. Gaiser (Gaiser, J., "Transverse isotropic Velocity Estimates From 
Slowness and Displacement Measurements", Ph.D. Thesis . 1989, University of Texas at Dallas, the 
teachings of which are incorporated by reference) did exactly this, using multi-offset (to perturb X 0 ) and 
multi-depth (to perturb Z 0 ) Vertical Seismic Profiles (VSP's). Gaiser used sources located on the surface 
and receivers located in vertical boreholes. 

In the present context, the luxury of a horizontal array of (fixed) sources and a vertical array of 
receivers is absent. Instead, the sonde 10 provides an array of sources and receivers (out to some 
maximum and minimum offsets) spaced along the borehole 9. The usual logging process notes that the 
source excites dispersive waves which are guided by the borehole walls 9 and which, for low frequencies, 
propagate with the medium's phase velocity in the direction of the borehole 9. Thus, by using low 
frequencies and analyzing the first breaks, one can determine (using standard borehole logging 
techniques) a single point on the anisotropic slowness-surface plot. 
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With this insight, the probiem is then one of how to make measurements over a range of 
propagation directions, instead of just the single direction given by the borehole. An important concept of 
this invention is that one can increase the dimensionality of measurements made by a borehole tool 10 by 
making use of seismic energy that has been reflected off of geological surfaces with a known orientation 
(geological strata, faults, the ocean bottom, or the surface of the ground). 

This is a new circumstance that becomes possible because use is made of a deviated borehole, 
drilled at an angle to the geology and/or to the surface. 

By contrast, in the textbook case of layer-cake geology, seismic isotropy and a vertical borehole, 
the reflections from the layers provide no additional anisotropy information, because the ray-paths from 
source to receiver are al! parallel. 

In the genera! case with a tilted borehole, arbitrary layer orientations, and general anisotropy, the 
extra information is complicated to make use of, requiring an anisotropic inversion problem to solve for the 
anisotropy. 

However, in the simpler case of polar anisotropy considered here, with a vertical symmetry axis 
and with vertical or horizontal reflectors, there is a simple geometrical interpretation of the reflected 
arrivals. Because the reflectors are a symmetry plane of the anisotropy, the incident and reflected angles 
at the reflector are equal. This allows us to use a standard geometrical artifice — image reflecting the 
source position through the plane of the reflector — to create an "image source" (Sj mage , of FIG. 5). 
Kinematically, the arrivals can then be processed as if they were direct arrivals from the image source 
Sj mage ,. In other words, as long as the borehole 9 is not entirely perpendicular to the reflector, one can 
convert the textbook one-dimensional borehole geometry into an effective two-dimensional borehole 
geometry, by using the mirror-image of the source in place of the "real" source. 

By this means, the borehole geometry has been converted into something like Gaiser's geometry. 
Then, by perturbing the receiver position, one can measure phase slowness along the direction of the 
borehole. And, by perturbing the source position, one can measure phase slowness along the direction of 
the image of the borehole. Just as in Gaiser's situation, one is directly finding points on the phase- 
slowness plot, but using a non-orthogonal coordinate system to do it (i.e., using the vector pointing along 
the borehole direction, and the vector pointing along the direction of the image of the borehole). This 
process becomes more uncertain as these two vectors are more parallel. For layer-cake geology, the 
optimum dip of the well is 45 degrees, in which case the borehole and the image of the borehole are 
perpendicular. 
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From the foregoing description, it will be observed that numerous variations, alternatives and 
modifications will be apparent to those skilled in the art. Accordingly, this description is to be construed as 
illustrative only and is for the purpose of teaching those skilled in the art the manner of carrying out the 
invention. Similarly, it should be understood that the disclosure, and claims which follow, are 
exemplifications of the principles of the invention. There is no intent to limit the invention to any particular 
set of anisotropy parameters, moveout equations, or mathematical algorithms. Almost any equation that 
couid be applied to surface geometry, could be applied to the environment of the present invention. With 
enough data, one might be able to get all 21 elastic constants of a general anisotropic medium. For 
example, if there were several fault planes nearby, with different strikes and dips, one could probe the 
medium with raypaths along a wide variety of azimuths and dips. That might be enough to determine the 
complete 3-dimensional local anisotropy. Moreover, equivalent computations may be substituted for those 
illustrated and described. For example, equation (11) of Tsvankin 1995 ( supra ) can be used in lieu of 
present Equation (8) when the reflector's normal coincides with the direction of symmetry and the moveout 
velocity is a function of the dip angle. Similarly, the more exact equations (10 through 17) of Thomsen 
1986 (supra ) can be used in lieu of present Equation (4). Also, the more general equations of Tsvankin 
and Thomsen 1994 ( supra ) and Tsvankin 1995 ( supra ) may be used to replace and augment the 
particular equations disclosed herein. Thus, it will be appreciated that various modifications, alternatives, 
variations, and changes may be made without departing from the spirit and scope of the invention as 
defined in the appended claims. It is, of course, intended to cover by the appended claims all such 
modifications involved within the scope of the claims. 
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CLAIMS 

We claim: 

1. A high resolution method for measuring seismic polar anisotropy in-situ in a region of 
interest characterized by geological formation elements that reflect sonic energy, comprising the steps of: 

(a) using a source to input sonic energy into the subsurface of the region of interest; 

(b) recording sonic energy reflected from the subsurface of the region of interest by 
using a plurality of sonic receivers that are located in a borehole that is located in the subsurface of the 
region of interest; said borehole having an axis that is at a known acute angle relative to at least one of the 
geological formation elements of said region; and 

(c) interpreting said recorded sonic energy of said receivers in terms of both the 
direct raypaths from said source to said receivers and the indirect reflected raypaths from said source to 
said receivers to obtain a measure of the seismic polar anisotropy parameters Vq.ii, and 6. 

2. The method of Claim 1 , wherein steps (a) and (b) are performed using a sonic energy 
source and a plurality of receivers that are located within said borehole. 

3. The method of Claim 2, wherein said source and said receivers are carried by a sonde 
that is adapted for axial movement within said borehole. 

4. The method of Claim 3, wherein said axis of said borehole is at an acute angle relative to 
the surface above said region; and wherein said geological formation elements comprise at least one of a 
bedding horizon, a fault, and a layer boundary having a generally horizontal orientation relative to said 
surface above said region. 

5. The method of Claim 3, wherein said geological formation elements comprise at least one 
of a dipping bedding horizon, a dipping fault, and a dipping layer boundary. 

6. The method of Claim 4, wherein the indirect raypaths of step (c) include reflections from 
geological formation elements comprising a plurality of iithologic horizons that are traversed by said 
borehole at an angle that is greater than zero degrees and less that ninety degrees, and that is preferably 
about 45 degrees; and wherein said axis of said borehole is deviated from the vertical. 

7. The method of Claim 6, wherein step (c) comprises the step of measuring the anisotropic 
phase slowness over a two-dimensional suite of directions, thereby characterizing the anisotropy of that 
part of the subsurface traversed by said raypaths. 
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8. The method of Claim 7, wherein said anisotropic phase slowness is measured over a two- 
dimensional suite of directions by at least performing the following steps: 

(a) perturbing the positions of said receivers, and 

(b) measuring phase slowness along said axis of said borehole. 

9. The method of Claim 7, wherein said anisotropic phase slowness is measured over a two- 
dimensional suite of directions by at least performing the following steps: 

(a) perturbing the position of said source, and 

(b) measuring phase slowness along the direction of an image of said borehole. 

10. The method of Claim 7, wherein step (c) is performed by finding points on a plot of phase- 
slowness, and using a non-orthogonal coordinate system to do so. 

11. The method of Claim 10, wherein said non-orthogona! coordinate system is defined by a 
vector pointing along said axis of said borehole, and by a vector pointing along the direction of an image of 
said borehole. 

12. The method of Claim 1 1 , wherein said borehole is inclined at about 45 degrees, such that 
said borehole and said image of said borehole are generally perpendicular to each other. 

1 3. The method of Claim 7, where in step (b) said recordings comprise waveforms excited by 
said source; and wherein step (c) comprises the steps of: 

(a) plotting said recorded waveforms as a function of time t and each source-receiver offset 

s; 

(b) identifying, on said plot of recorded waveforms, the arrival times of equal phase-points; 

and 

(c) measuring the slope ds/dt of a curve connecting said arrival times to determine the 
apparent phase velocity in the s-direction, whereby the inverse (dt/ds) of said slope is the corresponding 
phase slowness. 

14. The method of Claim 13, wherein said waveforms are direct reflected waveforms; and 
wherein said equal phase points are the peaks of corresponding arrivals. 

15. The method of Claim 6, wherein step (c) includes the step of computing 
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V p (9) = V 0 [l + 8 sin 2 e + ri sin 4 e] 



where V p (0) is the sonic P-wave velocity as a function of the angle of propagation 0 with respect to the 
5 symmetry axis of the medium defined by said lithologic horizons. 

16. The method of Claim 6, wherein said raypaths include: 

a. a direct borehole-parallel path, with arrival times that are a function of 

t, and 

10 where is the arrival time at the first receiver, X is the offset of a given receiver measured from the 
source and measured parallel to said borehole, and V p (9) is the sonic P-wave velocity as a function of the 
angle of propagation 0 with respect to the symmetry axis of the medium defined by said lithologic 
horizons. 

b. an indirect normal-incidence path, with arrival times that are a function of: 

x 

15 t 0 and — 

mo 

where t 0 is the time for normal-incidence reflection for the depth from said source to said reflectors 
comprising geological formation elements and V mo is the moveout velocity; and 

c. an indirect path of reflections from said lithologic horizons, with arrival times that 
are a function of: 

x 



20 1 0 , — — , and sin 0 



17. The method of Claim 16, the anisotropy parameter 5 and the vertical velocity V 0 are 
related to said moveout velocity by: 

V mo = V 0 Vl + 25 

mo 0 

18. The method of Claim 16, wherein the anisotropy parameter 5 and the vertical velocity V 0 
are reiated to said moveout velocity by: 
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V m „ = V 0 (1 + 5) 



19. The method of Claim 16, wherein said plurality of lithoiogic horizons are dipping; and 
wherein the anisotropy parameter 5 and the vertical velocity V 0 are related to said moveout velocity by: 



where 4» is the dip angle of the lithoiogic horizon traversed by said sonde in said borehole. 

20. Apparatus for measuring seismic poiar anisotropy in-situ in a region of the earth that is 
characterized by geological formation elements that reflect sonic energy and that comprise a plurality of 
lithoiogic horizons, comprising: 

(a) a housing that is adapted to travel within a borehole having at least one section that has 
an axis that is deviated from the vertical by a known acute angle and that is located at a known depth, said 
housing carrying at ieast one source of acoustic energy into the geological formation elements and a 
plurality of receivers for receiving acoustic energy from geological formation elements and said source; 

(b) a recording system, comprising elements that are located at the wellhead of said borehole 
and that are in communication with said source and said receivers, for receiving data therefrom; and 

(c) means for operating said source and said receivers and for processing data from said 
recording system in terms of both the direct raypaths from said source to said receivers and the indirect 
raypaths from said source, through the lithoiogic horizons that are traversed by said at least one section of 
said borehole, and to said receivers to obtain measures of at ieast the seismic polar anisotropy 
parameters V 0 , r\, and 5, said indirect ray-paths including reflections from geological formation elements 
that surround said borehole and at angles different from that of said axis of said at least one section of 
said borehole, said processing means comprising means for measuring the anisotropic phase slowness 
over a two-dimensional suite of directions. 

21. A method for measuring seismic polar anisotropy in-situ in a region of interest 
characterized by geological formation elements that reflect sonic energy, comprising the steps of: 



source being adapted for movement in a borehole that penetrates the region of interest, that is deviated 
from at least one of the vertical and said formation elements by a known non-zero acute angle; 



V 



m o 




(a) 



using a sonic energy source to input sonic energy into the region of interest, said 
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(b) recording sonic energy from said source and from said geological formation 
elements by using a plurality of axially spaced apart receivers that are adapted for movement in said 
borehole; and 

(c) interpreting the recorded signal in terms of both direct raypaths and indirect 
ray paths to measure at least the seismic polar anisotropy parameters V 0 , 5, and rj, said indirect ray- 
paths comprising a plurality of reflections from said geological formation elements that surround said 
borehole and reflections at angles different from that of the axis of said borehole, 

wherein the anisotropy parameter 5 and the vertical velocity V 0 are related to the moveout 

velocity V m0 by: 



v mo =v 0 VTT26. 

wherein the anisotropy parameters 8 and r\ and the vertical velocity V 0 are related to the P-wave 
velocity by: 

V p (0) = V o [l + 5 sin 2 9 + rj sin 4 0] 

15 where 9 is the angle of propagation with respect to the symmetry axis of the medium defined by said 
lithologic horizons; 

wherein the spacing X of said plurality of receivers located along the axis of the borehole, the 
P-wave velocity, and the arrival time ti at the first receiver are related to the arrival times t at the 
remaining receivers by: 

_ X 

20 1 ~ ^ + \T(0) ,and 

wherein said spacing X of said receivers, said moveout velocity, and the dip moveout time to are 
related to said arrival times t by: 

2/ x 2 x 2 2t n xsin0 

t 2 (x)= to + — r-+ °~ 

v } 0 v 2 V 

mo mo 

25 22. A high resolution method for measuring seismic polar anisotropy in-situ in a region of 

interest, comprising the steps of: 
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(a) inputting sonic energy into the subsurface by using a source located in a borehole 
that penetrates the region of interest and that is deviated from the vertical by a known non-zero acute 
angie; 

(b) recording sonic energy at a set of receivers located in said borehole; and 

5 (c) interpreting the received signal in terms of both direct raypaths and indirect 

raypaths, wherein said indirect raypaths include reflections from elements within the formations 
surrounding the borehole and at angles different from that of the axis of the borehole, said indirect 
raypaths including reflections from lithologic horizons. 

10 23. The method of Claim 22, wherein step (c) includes the steps of: measuring the seismic 

polar anisotropy parameters V 0 , 5, and r\; by computing: 

V p (9) = V 0 [l + 5 sin 2 0 + ri sin 4 9] 

where V (6) is sonic P-wave velocity as a function of the angle of propagation 0 with respect to the 
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p 

symmetry axis of the medium. 

24. The method of Claim 22, wherein said source and receiver are in a single borehole; and 
said raypaths include: 

a. a direct borehole-parallel path, with arrival times given by 

x 

t = t, + 7—7" 

' v,(e) 

20 where t , is the time to the first receiver, and offset x is measured from that receiver; and 

b. an indirect normai-incidence path, with arrival times given by 

t 2 ( x ) =t o+^ , 



mo 



where t 0 is the time for normal-incidence reflection for the depth from the source to the reflector; and 

c. an indirect path of reflections from horizontal interfaces, with arrival times given 

25 by: 

mo mo 
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where offset at is the offset of a given receiver, measured from the source, parallel to the borehole, and 
V mo , the moveout velocity of such reflected arrival, is related to the desired anisotropy parameters by: 

v mo =v 0 VTT25. 



25. The method of Claim 24, wherein V 0 Vl + 25 is approximated by V 0 [l + 5] . 

26. The method of Claim 22, wherein step (c) comprises the step of measuring the 
anisotropic phase slowness over a two dimensional suite of directions to characterize the anisotropy of the 
medium traversed by the rays, said step of measuring the anisotropic phase slowness including the steps 
of: 

(a) perturbing the position of at least one of said receivers and said source; and 

(b) measuring phase slowness along said axis of said borehole. 

27. The method of Claim 26, wherein wherein step (c) is performed by finding points on a plot 
of phase-slowness, and using a non-orthogonal coordinate system defined by a vector pointing aiong said 
axis of said borehole, and by a vector pointing along the direction of an image of said borehole. 

28. The method of Claim 22, where in step (b) said recordings comprise waveforms excited 
by said source; and wherein step (c) comprises the steps of: 

(a) plotting said recorded waveforms as a function of time t and each source-receiver offset 

s; 

(b) identifying, on said plot of recorded waveforms, the arrival times of equal phase-points; 

and 

(c) measuring the slope dt/ds of a curve connecting said arrival times to determine the 
apparent phase velocity in the s-direction, whereby the inverse (ds/dt) of said slope is the corresponding 
phase slowness. 

29. The method of Claim 28, wherein said waveforms are direct reflected waveforms; and 
wherein said equal phase points are the peaks of corresponding arrivals. 

30. A seismic method, comprising the steps of: 

(a) locating a borehole tool at a non-zero acute angle to the horizon and measuring said 
acute angle of inclination of said borehole tool, said borehole tool carrying plurality of receivers that are 
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located along the axis of said borehole and that are spaced apart from one another, said borehole tool 
carrying a seismic source that is spaced apart from said receivers ; 

(b) measuring the spacing of each receiver from said source along said axis; 

(c) transmitting seismic energy from a source carried by said borehole tool, and measuring 
the arrival time at a first receiver that is closest to said source and the arrival times at the remaining 
receivers; 

(d) obtaining a mathematical relationship between the anisotropy parameter 6* , vertical 
velocity and moveout velocity; 

(e) obtaining a mathematical relationship between the anisotropy parameters 8 and r\, said 
vertical velocity, said angle of inclination of a borehole tool and P-wave velocity; 

(f) obtaining a mathematical relationship between said spacing of said receivers along the 
axis of the borehole, P-wave velocity, said arrival time at a first receiver and said arrival times at said 
remaining receivers; 

(g) obtaining a mathematical relationship between said spacing of said receivers, said 
moveout velocity, said angle of inclination, dip moveout time, and said arrival times; and 

(h) using the measurements of steps (a), (b), (c) and the mathematical relationships of steps 
(d), (e), (f) and (g) to determine measures of at least said seismic anisotropy parameters 6 and n- 

31. In seismic exploration of a region of the earth containing at least one of a bedding horizon, 
a fault, and a layer boundary of a known orientation, wherein a borehole tool is used that carries a plurality 
of receivers that are located at positions x along the longitudinal axis of the tool and that are spaced apart 
from one another, and that carries a seismic source that is spaced apart from the receivers, a process 
comprising the steps of: 

(a) locating the borehole tool in a borehole that is inclined to at least one of the orientation of 
the bedding horizon, the fault and the layer boundary; 

(b) recording waveforms from direct raypaths from the source to the receivers and from the 
indirect raypaths from the source through at least one of a bedding horizon, a fault, and a layer boundary 
of a known orientation; 

(c) obtaining a mathematical relationship between the anisotropy parameters 6 and n, the 
vertical velocity V 0 , the angle of propagation 8 with respect to the symmetry axis of the region traversed 

by the borehole, and the P-wave velocity V p ; 

(d) obtaining a mathematical relationship between the spacing of the receivers along the axis 
of the borehole, said angle of inclination of the borehole, dip moveout time, and said arrival times; and 
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(e) using the recorded waveforms of step (b) and the relationships of steps (c) and (d) to 
obtain the anisotropy parameters V 0 , 5, and r\. 

32. The seismic method of claim 31, where V mo is a function of V 0 and 8. 



33. The seismic method of claim 32, where V mo — V 0 Vl + 28 . 

34. The seismic method of claim 31, where V p is a function of 5, rj, and 0. 

2 

10 35. The seismic method of ciaim 34, where V p is a function of V 0 , 8, rj, and Sin 0. 

36. The seismic method of claim 35, where 



V p = V 0 [l + 5 sin 2 e + rj sin 4 0] 



15 



37. The seismic method of claim 31, where said arrival times are a function of x, V mo , 0, 
and the vertical travel time tr> 

38. The seismic method of claim 37, where said arrival times are a function of X, V mo , sinG, 
20 and the vertical travel time Xq. 

39. The seismic method of claim 38, where! 

2 2 x 2 2t n xsin0 

t 2 = t 2 + T-+ ~ 

0 vl V, 



mo mo 



25 40. The seismic method of Claim 31, wherein said at least one of a bedding horizon, a fault, 

and a layer boundary of a known orientation are dipping; and wherein the anisotropy parameter 8 and the 
vertical velocity V 0 are related to said moveout velocity by: 



21 



WO 98/59264 



PCT/US98/12866 



m o 



cos§ V V 0 ((()) dQ 



where <j> is the dip angle of said at least one of a bedding horizon, a fault, and a layer boundary traversed 
by said source and said receiver. 
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